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ABSTRACT 



We perform photometric analyses for a bright early-type galaxy (ETG) sam- 
ple with 2949 galaxies (Mr < —22.5 mag) in the redshift range of 0.05 to 0.15, 
drawn from the SDSS DR7 with morphological classification from Galaxy Zoo 1. 
We measure the Petrosian and isophotal magnitudes, as well as the correspond- 
ing half-light radius for each galaxy. We find that our Petrosian magnitudes, 
and isophotal magnitudes to 25 mag/arcsec^ and 1% of the sky brightness are 
on average 0.16 mag, 0.2 mag, and 0.26 mag brighter than the SDSS Petrosian 
values, respectively. In the first case the underestimations are caused by overesti- 
mations in the sky background by the SDSS PHOTO algorithm, while the latter 
two are also due to deeper photometry. We find the overestimations are more 
strongly influenced by galaxy sizes than by galaxy luminosities. Similarly, the 
typical half-light radii {r^o) measured by the SDSS algorithm are smaller than our 
measurements. As a result, the bright-end of the r— band luminosity function is 
found to decline more slowly than previous works. Our measured luminosity den- 
sities at the bright end are more than one order of magnitude higher than those 
of Blanton et al. (2003), and the stellar mass densities at ~ 5 x lO^^M© and 
M* ~ are a few tenths and a factor of few higher than those of Bernardi 

et al. (2010). These results may significantly alleviate the tension in the assem- 
bly of massive galaxies between observations and predictions of the hierarchical 
structure formation model. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: luminosity function, 
mass function - galaxies: photometry 
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1. INTRODUCTION 

The properties of nearby early type galaxies (ETGs), especially the brightest ones, 
offer important clues to understanding the cosmic assembly history of massive galaxies. In 
terms of their morphologies, colors, stellar population content and scaling relations, ETGs 
appear to be relatively simple systems compared with spirals and other galaxies. There 
are, however, a lot of renewed interests in their dynamical properties, particularly from the 
recent integral field unit surveys (e.g. Emsellem et al. 2007; Brough et al. 2010; Cappellari 
et al. 2012). There is also a debate on their stellar mass assembly processes (Renzini 2006; 
Scarlata et al. 2007), and so the understanding of ETGs remains a particularly interesting 
issue. 

In the past two decades, the concept of "downsizing" (Cowie et al. 1996; Gavazzi 
& Scodeggio 1996; Fontanot et al. 2009) for galaxy formation has been widely discussed. 
According to this scenario, the epoch of star formation in ETGs depends on the galaxy 
mass, namely, the most massive ETGs formed their stars on a shorter time scale and at 
earlier times. Furthermore, large ground-based and space-based imaging and spectroscopic 
surveys at low and high redshifts found that the mass function shows a weak evolution for 
massive ETGs since redshift ~ 1 (e.g. Cimatti et al. 2006; Scarlata et al. 2007; Cool et 
al. 2008; Cimatti 2009). In contrast, theory predicts that the typical stellar mass of most 
massive ETGs increases by a factor of 2-4 since redshift one (see De Lucia & Blaizot 2007; 
Tonini et al. 2012). 

The size evolution for ETGs has also been the topic of many recent studies. It has 
been shown that the massive and quiescent ETGs with stellar mass > lO^^M© at high 
redshift are more compact with effective radii a factor of ~ 3-5 smaller than the present-day 
ETGs with similar stellar masses, but the stellar masses increase by a factor of ~ 2 since 
redshift 2 (e.g. Daddi et al. 2005; Trujillo et al. 2006; van Dokkum et al. 2008). A popular 
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explanation for the evolution of massive ETGs is minor dry mergers, although additional 
physical mechanisms may be required (e.g. Cimatti et al. 2008; Bezanson et al. 2009; 
van Dokkum et al. 2010; Bruce et al. 2012 and McLure et al. 2012). It should be noted 
that such studies have not focused on the evolution of ETGs at the high mass end with 
M* > 1O^^M0. As pointed out by Naab (2012), the most massive ETGs or their progenitors 
start forming their stars at redshift ~6 or even earlier and these high redshift galaxies are 
just the cores of their local counterparts. The assembly at later times may be inside out 
by accretion of stellar mass in the outskirts of these galaxies. However, the mass assembly 
history of most massive ETGs is not yet completely clear: the current ACDM model cannot 
reproduce the stellar age, metallicity and color evolution simultaneously (Tonini et al. 2012; 
De Lucia & Borgani 2012). 

To understand the mass and size evolution, accurate photometry for massive ETGs is 
an essential first step. Bright galaxies tend to be located in crowded environments or have 
extended stellar halos. Inadequate masking of their neighbors or extended halos can result 
in an overestimate of the sky background, leading to an underestimate of their luminosity 
and size. Indeed a number of investigators have noted that the SDSS photometric reduction 
systematically underestimates the luminosities and half-light radii for bright ETGs (see 
Aihara et al. 2011 for the SDSS DR8 release and DR7 documentation). For the SDSS 
DR8, a more sophisticated sky background subtraction algorithm has been adopted and 
significant improvements were achieved, but some problems still persist (see §4.1). On 
the other hand, the surface brightness distribution of most massive ETGs cannot be fit 
well by a simple model such as the commonly used Sersic (1963) law, which can also lead 
to an underestimate of the luminosities. Graham et al. (2005) carefully investigated the 
deficiency of Petrosian (1976) system to quantify galaxy luminosity and size. They find 
that the Petrosian magnitude and size commonly used by studies based on SDSS strongly 
depend on the surface brightness profile. For example, the flux deficiency can be as large 
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as 0.5 mag for a galaxy with r^/^ surface brightness profile. To avoid underestimating 
the luminosities of the brightest cluster galaxies (BCGs), von der Linden et al. (2007) 
used isophotal photometry to 23 mag/arcsec^ in the r— band, rather than the Petrosian 
magnitude. Therefore, aperture photometry offers a viable alternative to describe the 
photometric properties of ETGs. 

Underestimating the luminosities of luminous ETGs will lead to an underestimate of 
the stellar mass density at the bright end. Using their their own sky background subtraction 
algorithm (cmodel) and stellar mass estimation, Bernardi et al. (2010) found a higher 
number of very massive galaxies than previous works. The excess can be up to a factor of 
~ 10 when the stellar mass M* is larger than 5 x IQi^^Mq, which is highly significant. 

In this work, we perform accurate photometry for a complete sample of nearby bright 
early-type galaxies with My. < —22.5 mag and > IO^^Mq. The morphologies of these 
galaxies are taken from the Galaxy Zoo project which classified nearly 900,000 SDSS 
galaxies (Lintott et al. 2011). We perform our own sky background subtraction, which 
leads to much more accurate photometry of ETGs. We then compare the luminosities and 
sizes measured by different methods and investigate the luminosity function as well as the 
stellar mass density for the bright ETGs. 

The structure of this paper is as follows. In §2 and §3 we describe the sample selection 
and data reduction for our bright ETGs. In §4, we compare our measured luminosities and 
sizes of bright ETGs with the SDSS Petrosian magnitudes and sizes, and present the results 
of the r— band luminosity function and the stellar mass density. We finish the paper with 
a summary in §5. In this paper, we adopt a Hubble constant of Hq = 70kms~^Mpc~^, 
matter density parameter Vt^ = 0.3 and cosmological constant Vt\ = 0.7. 
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2. SAMPLE 

Our bright early type galaxy (ETG) sample is drawn from the morphological catalogue 
of Galaxy Zoo 1 (Lintott et al. 2011). The Galaxy Zoo 1 project performed visual 
morphological classifications for nearly 900,000 galaxies from the Sloan Digital Sky Survey 
(SDSS, York et al. 2000), in which 667,945 galaxies are from the Main Galaxy Sample 
of SDSS (MGS, Strauss et al. 2002). The MGS includes galaxies with r— band Petrosian 
magnitude r < 17.77 and r — band Petrosian half-light surface brightness fi^Q < 24.5 
mag/arcsec^. All 667,945 galaxies from the MGS have spectroscopic redshift in the range 
of 0.001 < z < 0.25 and u, g, r, i, z band photometry based on SDSS DR7 (Lintott et al. 
2011). 

Given that the SDSS spectroscopy survey is incomplete for bright galaxies with redshift 
less than 0.05 (Stoughton et al. 2002; Strauss et al. 2002; Schawinski et al. 2007; Kaviraj 
et al. 2007) and reliable photometric analyses with high signal to noise (S/N) ratios can 
only be performed for galaxies with r < 16 mag (Fukugita et al. 2007), we constrain our 
bright ETG sample to z > 0.05 and r < 16 mag. Since we are only concerned with the 
photometric properties of the luminous ETGs, we further restrict our sample to bright 
ETGs with Mr < —22.5 mag. The Petrosian absolute magnitudes, M^, are calculated using 
the equation Mr = mr — 5 log(i5L/10 pc) — A — k, where Dl is the luminosity distance, A is 
the Galactic extinction obtained from the photometric catalogue of SDSS DR7 and k is the 
k-correction derived using the IDL KCORRECT algorithm of Blanton et al. (2007). Finally 
there are 7930 bright ETGs in the redshift range of 0.05 < z < 0.15 and Mr < —22.5 mag. 

We further divide this sample into three volume-limited subsamples with 2303 ETGs in 
the redshift ranges of 0.05 < ^ < 0.1 for Mr < -22.5 mag, 538 ETGs with 0.1 < z < 0.125 
for Mr < -23 mag and 108 ETGs with 0.125 < z < 0.15 for Mr < -23.5 mag, as shown 
in Fig. [H Then we perform < V/Vmax > test (Schmidt 1968) for these subsamples. The 
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average value < V^/Vmax > are 0.51, 0.50 and 0.51 for the three subsamples, respectively, 
which indicates that all the three subsamples are spatially homogeneous. The photometric 
analysis in this work is based on these three subsamples. In total, there are 2949 ETGs 
brighter than —22.5 mag, among which 1053 ETGs are brighter than —23 mag. Fig. [2] 
shows the spatial distribution of our sample galaxies on the sky in Galactic coordinates: 
the total area is 9055 degree^, corresponding to 22% of the whole sky. Fig. [3] presents the 
histograms showing the distributions of apparent magnitude and spectroscopic redshift 
for 2949 early- type galaxies. The median Petrosian magnitude is 15.37 mag and median 
redshift is 0.087, respectively. 

In order to verify the early-type morphology of our sample galaxies, we visually 
inspected the SDSS r— band images for all 2949 galaxies. In addition, we also checked other 
commonly used classification criteria for early-type galaxies. All our sample galaxies have 
g — r > 0.7 (e.g. Blanton et al. 2003; Shen et al. 2003), and most (94%) have concentration 
index Cr = rgo/r^o larger than 2.86 (Bernardi et al. 2010), where rgo and rso are the radii 
containing 90% and 50% of Petrosian flux. Therefore, our selected sample galaxies are all 
ETGs. 

3. DATA REDUCTION 

3.1. Estimation of the Sky Background 

The corrected frame fpC-images in the r— band for our sample galaxies are directly 
obtained from the SDSS DR7 Data Archive Server. The images have been pre-processed 
by the SDSS photometric pipeline (PHOTO), which includes bias subtraction, fiat-fielding 
and bad pixels correction (cosmic rays removal, bad columns and bleed trails). In order to 
obtain the SDSS photometric flux calibration information during observations, such as the 
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photometric zeropoint a , the first-order extinction coefficient k and the airmass X, we also 
downloaded the calibrated field statistic file, named as tsField from the archive. 

As mentioned in §1, the SDSS photometric reduction systematically underestimates 
the luminosities and half-light radii for bright ETGs, which is mainly caused by inadequate 
masking of their neighbors or extended stellar halos (Aihara et al. 2011), leading to an 
overestimate of the sky background. In this work we perform the masking more carefully 
and estimate the sky background model following Liu et al. (2008), which has been 
successfully used to measure the luminosities and half-light radii for the brightest cluster 
galaxies (BCGs) that are located in crowded fields and usually have extended faint stellar 
halos. 

In the following, we outline our sky background subtraction approach. In order to 
obtain the sky background, we first masked out all objects detected by SExtractor (Bertin 
& Arounts 1996) in the corrected frames with 2048 x 1489 pixels (13'. 5 x 9'. 8). We then 
carefully checked each frame by eye to make sure that the wings of bright stars or the 
faint stellar halos of galaxies have been properly masked. As expected, we find that the 
automatic algorithm of SExtractor does not work well for about 15% bright ETGs, which 
mostly include ETGs residing in crowded fields, having extended stellar halos or are close 
to foreground bright stars. For the 2949 ETGs brighter than —22.5 mag, the percentages of 
these objects are about 13%, 7%, 4%, respectively. For the 1053 ETGs brighter than —23 
mag, the percentages are even higher, about 18%, 10%, 5%, respectively. For these ETGs, 
we modify the mask images manually. For comparison. Fig. H] shows three cases (top, middle 
and bottom rows) between the masked images generated by the automatic algorithm of 
SExtractor (middle panels) and those by hand (right panels). The left panels of Fig. SJshow 
the original true color images for the example galaxies. The top row shows a galaxy in a 
crowded field. From the top middle panel, it can be seen that SExtractor cannot separate 
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the target galaxy from its neighbors well. Therefore we manually flag the nearby objects 
with circles that are large enough to cover the objects completely, as shown in the top right 
panel. The second row shows a bright ETG with an extended stellar halo. From the middle 
panel, it is clear that the flagged area generated by SExtractor is not big enough to cover 
the whole stellar halo. Hence, the stellar halo would be considered as the background and 
subtracted from the galaxy itself, leading to underestimates of the luminosity and half-light 
radius of the galaxy. Following the shape of the target ETG, we use an ellipse to mask the 
entire stellar envelope, as shown in the middle right panel. SExtractor also cannot handle 
well galaxies surrounded by bright foreground stars, especially those with diffractive spikes 
(bottom panels). Therefore, we use long rectangles to mask the star spikes (bottom right). 
The masked images generated manually can provide not only more accurate sky background 
images but also good masks for the surface photometry (see §4) . 

Next we smooth the sky background only image with a median filter of 51 x 51 pixels. 
After the median filtering is performed, masked regions smaller than the median box filter 
are replaced with the surrounding sky background, whereas part of the masked regions 
larger than the median box filter remain flagged. Then a two-dimensional second-order 
Legendre polynomials is fitted to the filtered sky background image using IRAF/IMSURFIT 
task. Finally, a Gaussian filter with a standard deviation a = 9 pixels is convolved with the 
fitted frame to obtain the final sky background image. The sky background map is typically 
tilted with a spatial variation of ~ 1 — 2 ADU across the whole frame (Liu et al. 2008). We 
subtract this sky background model from the initial fpC-image to obtain the sky-corrected 
frame. The blank regions in the sky-subtracted frame follow Gaussian distributions with 
means close to zero and standard deviations of several ADU, which is consistent with Liu et 
al. (2008). After the sky subtraction, we trim both the sky-subtracted frame and the mask 
image to 501 x 501 pixels centered on the galaxy of interests. The trimmed mask image 
with the target galaxy un-flagged will be used to probe the masked regions in the isophote 
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fitting. 

3.2. Isophotal Photometry 

After the sky background subtraction, we perform surface photometry for our bright 
ETGs in order to estimate the luminosities and sizes of sample galaxies. As is well known, 
the Petrosian magnitude and Petrosian size are most commonly used to describe the galaxy 
flux and half-light radius for the analyses based on SDSS database, because they do not 
depend on the model fitting to galaxies. The Petrosian radius Vp is defined as the radius r 
at which the ratio of the local surface brightness averaged over an annulus between 0.8r and 
1.25r to the mean surface brightness within r equals to 0.2 and the Petrosian magnitude 
is the integrated flux within 2rp (Petrosian 1976). However, the Petrosian magnitude 
misses the light outside 2rp (Petrosian aperture) because of its dependence on the surface 
brightness profile of galaxies (Graham et al. 2005), which leads to underestimates of the 
fluxes and sizes for galaxies with extended stellar halos. In this work, we will measure 
not only the Petrosian magnitude and Petrosian half-light radius, but also the isophotal 
magnitude and half-light radius to deeper isophote limits at 25 mag/arcsec^ and 1% of sky 
brightness, corresponding to ~ 26 mag/arcsec^ in most cases (Bernardi et al. 2007). 

We perform the surface photometry analysis following Wu et al. (2005) and Liu et al. 
(2008). The procedures are briefly described below. First, we use ISOPHOTE/ELLIPSE task 
in IRAF to flt each of the trimmed sky-subtracted images excluding the masked regions in 
the fitting. The surface brightness of the target galaxy is fitted by a series of elliptical annuli 
in a logarithmic step of 0.1 along the semi-major axis. The annuli chosen in the outer parts 
of image is larger, which can suppress the shot noise in the outer regions where the signal 
to noise ratio (S/N) is much lower. The output of ISOPHOTE/ELLIPSE is the mean intensity 
in each isophote annulus. Then we integrate the surface brightness profile to isophotal 
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limits of 25 mag/arcsec^ or 1% of sky brightness to obtain the apparent magnitudes and 
half-hght radius r^Q. We also measure the Petrosian magnitudes and rso based on our sky 
background subtracted images. In our analysis, the equivalent radius \/ab of an ellipse is 
used for all the radial profiles, where a and b are the semi-major and semi-minor radii of 
the ellipse. The cosmological dimming is also taken into account for the surface brightness 
profiles. The observational errors in the surface brightness profile include random errors 
(e.g., the shot noise of the object and sky background, readout noise and noise contributed 
by data reduction) and the error from sky background subtraction. 

4. RESULTS 

4.1. Luminosities of bright ETGs in the local universe 

The SDSS database provides the largest galaxy sample with both photometric and 
spectroscopic information in the local universe. However, the underestimation of the 
luminosities and sizes for the bright ETGs prevents us from correctly understanding their 
properties and assembly history. 

For brevity in later discussions, we define 

Amp = mp - mp^sdss7, 

Am25 = - mp^sdss7, (1) 

Ami% = mi% - mp,sdss7, 
where mp, m25 and mi% are our measured Petrosian magnitude, isophotal magnitudes with 
surface brightness measured to 25 mag/arcsec^ and 1% of sky brightness, and mp sdss7 is the 
Petrosian magnitude from the SDSS DR7 pipeline in the r— band. 

The top and middle panels in Fig. [5] show Amp (defined eq. [1]) as a function of 
the apparent and absolute magnitudes, while the bottom panel shows the histogram of 
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Amp. Clearly for more luminous ETGs, the luminosity difference between the SDSS and 
our measurements is larger. The mean and median values of the luminosity differences 
are 0.16 mag and 0.14 mag, respectively. It shows that the algorithm used in the SDSS 
DR7 has systematically overestimated the sky background for bright ETGs, leading to 
underestimates in the luminosities of bright ETGs, especially for the brightest ones. 

Addressing the same issue, in the SDSS DR8, Aihara et. al. (2011) have re-processed 
all SDSS imaging data using a more sophisticated sky background subtraction algorithm 
and obtained significant improvement. The left and right panels of Fig. [Hlare histograms 
of luminosity difference between the SDSS DR8 and SDSS DR7, and between our 
measurements and SDSS DR8, respectively. We can see that the median values of 
luminosity difference for these two cases are —0.05 mag and —0.08 mag, and the mean 
values are about —0.04 mag and —0.12 mag, respectively. The new algorithm on sky 
background subtraction used in the SDSS DR8 has indeed improved the photometric 
measurement (the mean shift is about 0.04 mag). However, compared to our study it 
still underestimates the luminosities by about 0.12 mag on average for bright ETGs. We 
carefully checked each image with very large luminosity differences between the SDSS DR8 
and our measurements. We find that most of these galaxies reside in crowded fields or have 
very extended faint stellar halos, and so the SDSS DR8 treatment of the sky background 
issue is incomplete. 

Liu et al. (2008) found that the isophotal magnitudes measured to the surface 
brightness of 25 mag/arcsec^ are generally larger than the Petrosian values for BCGs. We 
thus also compare the Petrosian and isophotal magnitudes for our sample ETGs. The top 
and middle panels of Fig. [7] shows Am25 (defined in eq. [1]), as a function of and M25, 
respectively. It is obvious that Am25 increases with the apparent and absolute magnitudes 
of ETGs. The mean and median differences are 0.20 mag and 0.17 mag, respectively, which 
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are larger than those for Arrip (0.16 and 0.14 mag, respectively). 

Bernardi et al. (2007) investigated the galaxy luminosity by integrating the best-fit 
model of galaxy surface brightness profile to 1% of sky brightness. Fig. [8] shows Ami% 
(defined in eq. [1]) as a function of mi% and Mi%, respectively. The same trend is seen: 
Ami% increases as ETGs become more luminous. The mean and median values of Ami% 
are 0.26 mag and 0.23 mag, respectively. Note that the photometric limit of 1% sky is ~ 26 
mag/arcsec^, almost one magnitude deeper than 25 mag/arcsec^. 

We further examine the percentage of ETGs with large luminosity difference between 
the SDSS Petrosian and our measured magnitudes. For the 2949 ETGs brighter than —22.5 
mag, there are 6%, 12% and 22% ETGs with Amp, Am25 and Ami% larger than 0.3 mag. 
For the 1053 ETGs brighter than —23 mag, the fractions are even higher, 11%, 19% and 
31%, respectively. For more luminous ETGs, the underestimation by the SDSS algorithm 
is more severe. In addition, the deeper the photometric measurements, the larger the 
underestimations. 

We point out that to the photometric limit of ~ 26 mag/arcsec^ the light still belongs 
to galaxies as shown by Tal & van Dokkum (2011). They stacked more than 42000 SDSS 
images of LRGs (Luminous Red Galaxies), reaching a depth of ~30 mag/arcsec^, and 
found that the stellar light out to 100 kpc is physically associated with galaxies, instead 
of inter-cluster or inter-group light. Our photometric measurement to 1% sky brightness 
reaches at most 100 kpc (50 kpc on average, see §4.2), and so our measured isophotal light 
is from ETGs themselves. 
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4.2. Sizes of bright ETGs in the local universe 

Galaxy size is one of the most important parameters of galaxy properties. A reliable 
determination of the sizes of bright ETGs in the local universe provides the basic calibration 
for investigating the size evolution, an area of active study in recent years (e.g. Szomoru et 
al. 2012; McLure et al. 2012; Trujillo 2012). The aforementioned works compare galaxy 
sizes at high redshift to those of the local universe by Shen et al. (2003) based on the SDSS 
PHOTO algorithm, who found a power-law relation between the galaxy luminosity and 
size (Shen et al. 2003). However, if the luminosities of galaxies have been underestimated, 
galaxy sizes may have been underestimated too. In this section, we will discuss the 
size difference between the SDSS and our measurements based on different luminosity 
estimations. For convenience, we define r5o,p, ^"50,25 and rso, 1% as our Petrosian half-light 
radius, isophotal half-light radii to 25 mag/arcsec^ and 1% of the sky brightness. These 
values will be compared with rso^sdss?; the Petrosian half-light radius from the SDSS DR7 
pipeline in the r— band. 

Fig. [9] present the histograms of the half-light radii distributions measured by different 
photometric methods. The median values of rso^sdss?, ''^50, ''^50,25 and ^50,1% are ~ 8.6, 10.4, 
11 and 12 kpc, respectively. We can clearly see from Fig. [9] that the largest r5o,sdss7 value is 
smaller than 20 kpc, while our measured r5o,p, tso, 25 and ^50,1% can be as large as 30 kpc 
and a large fraction (~ 27%) of brightest ETGs have sizes larger than 15 kpc. Therefore, 
the SDSS algorithm has significantly underestimated the real sizes of bright ETGs. 



We fit the size distribution with a log-normal distribution. 



/(r,r,ainr-) = -7=^ exp 



In^ (r/r) 
2<tL 



'A (2) 

r 



which is characterized by the median r and the dispersion o"inr- We find the best-fit 
values for rso^sdss?, '"so.p! "'^50,25, ^"50,1% distributions are f = 8.58, 10.28, 10.80, 11.82 kpc and 
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o"inr = 0.24,0.32,0.31,0.32, which are consistent with the median values mentioned above. 

Fig. ITU] shows the luminosity differences. Amp, Am25 and Ami%, as a function of r^Q p, 
1^50,25 and r5o,i%, respectively. We can see from Fig. [10] that the luminosity difference is 
well correlated with the size of ETGs measured by different approaches. It is interesting 
that the correlations between Am and r^Q shown in Fig. [TU] are tighter than the correlation 
between Am and luminosity as shown in Figs. [5], [7| and [H The statistical scatters in the 
correlations between Am and rso are much smaller than those in the correlations between 
Am and luminosities. This leads to a conclusion that the luminosity underestimation 
depends more strongly on the size of ETGs than on their luminosity. 

We further examine the images of the brightest ETGs (Mj. < —23 mag) with differences 
between the ri% and the SDSS DR7 measurement larger than 10 kpc and Ami% > 0.4 mag. 
In total, there are 102 such ETGs; 56% ETGs are in crowded field, 26% have extended 
stellar halos and the remaining 18% are contaminated by bright nearby stars. Fig. [TT] shows 
typical examples of ETGs located in crowded fields (top panel) and those contaminated 
by nearby bright stars (bottom panel). Fig. [T^ represents the ETGs with extended stellar 
halos (top panel) and corresponding surface brightness profiles and the residuals from 
Sersic models (bottom panel). It is obvious from Figs. [TT] and [12] that any automatic 
algorithm cannot do the sky background subtraction well for these ETGs, leading to a 
serious underestimate of the luminosity and size. 

Fig. [I3] shows the size and luminosity relations measured by different methods. In 
each panel, we give the best power-law fit slope for the correlations. The top left and 
right panels are for Petrosian half-light radius with Petrosian absolute magnitude obtained 
from the SDSS DR7 catalog and measured by us, respectively. The slope in the top right 
panel {a = 0.90 ± 0.03) is steeper than that in the top left panel {a = 0.76 ± 0.02). It 
indicates that the sky background subtraction for bright ETGs significantly influences the 
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size-luminosity relation. The bottom left and right panels of Fig. [13] show the correlations 
between the size and isophotal absolute magnitude to 25 mag/arcsec^ and to 1% of sky 
brightness, respectively. We can see from the bottom two panels of Fig. [12] that the slope for 
the correlation between logr5o,i% and Mi% {a = 0.92 ± 0.03) is somewhat steeper than the 
correlation between logr5o,25 and M25 (a = 0.87 ± 0.02). These results are consistent with 
the size-luminosity relation trend for the brightest cluster galaxies (BCGs) obtained by Liu 
et al. (2008) who found that the power-law slope becomes steeper when the measurement 
goes deeper (see also Bernardi et al. 2007). 

Our derived slope on the logrso and Mj. relation for ETGs is much larger than the 
value 0.67 found by Shen at al. (2003) that is widely adopted by recent works on the size 
evolution for bright ETGs. If we use our measured sizes of ETGs, the size evolution since 
redshift 2 will be somewhat larger. However, due to the surface brightness dimming, it may 
be difficult to perform photometry down to 26 mag/arcsec^ in the rest-frame r— band for 
ETGs at redshift 2. As a result, we should be more cautious in discussing the size evolution 
by explicitly taking into account the survey surface brightness limit. 



4.3. The bright end of the r— band LF and stellar mass density 

A basic way to investigate galaxy properties and their evolution is by studying the 
luminosity function (LF). There are already many LF studies using different samples and 
approaches at different redshifts (see the recent review paper by Johnston 2011). Given 
that the luminosities of the brightest galaxies (—23.5 mag) have been underestimated (by 
~10% to 40%), it is worth revisiting the LF at the bright end in the local universe. 

In this work, we construct the galaxy luminosity function at the bright end, utilizing 
the non-parametric 1/Knax method (Schmidt 1968; Felten 1976; Eales 1993). Our bright 
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ETG sample includes 2949 galaxies consisting of three subsamples as described in §2, in 
which we have already discussed the homogeneity for this sample. Briefly, the 1/Vmax is 
the inverse of the maximum volume, to which the galaxy could have been detected. The 
LF is obtained by integrating 1/Knax in different luminosity bins for the whole sample of 
galaxies. Given that our ETGs contain three subsamples, we calculate the LF in three 
volumes separately, and then average them to obtain the final LF. In addition, SDSS fiber 
collisions lead to ~ 7% incompleteness for the spectroscopic sample (Bernardi et al. 2010), 
we multiply the 1/Knax counts by a factor of 1/0.93 to obtain the final LF. 

The top panel of Fig. [H] shows our r— band luminosity function at the bright end. The 
green triangles, red circles and blue solid squares represent the luminosity function for our 
measured Petrosian magnitude, isophotal magnitudes measured to the surface brightness 
of 25 mag/arcsec^ and 1% of the sky brightness, respectively. For comparison, we also 
plot the luminosity functions from Blanton et al. (2003) and Bernardi et al. (2010) in 
the top panel of Fig. [TH Blanton et al. (2003) used a sample of 147,986 galaxies (from 
SDSS EDR) and the maximum likelihood method to calculate the luminosity function at 
z = 0.1. Their sample is much larger, but includes all morphological types. However, at the 
bright end (e.g. the galaxies brighter than —22.5 mag), their galaxies are almost exclusively 
ETGs. The Petrosian r— band luminosity of galaxies in Blanton et al. (2003) are directly 
obtained from the SDSS catalogue based on the SDSS PHOTO algorithm. Bernardi et al. 
(2010), on the other hand, used an ETG sample of galaxies selected from ~ 250000 SDSS 
galaxies with 14.5 < mp^sdss < 17.5 using the concentration index Cr > 2.86, which is a 
conservative way to select ETGs from SDSS. The luminosities of ETGs in Bernardi et al. 
(2010) are calculated using the cmodel in the SDSS pipehne with their own sky background 
subtraction method. 

It is clear from the top panel of Fig. [TUthat at the bright end (Mj. < —22.5 mag), the 
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slope of the Bernardi et al. (2010) 's LF is shallower than that of Blanton et al. (2003), 
implying that the luminosity estimation method of Bernardi et al. (2010) is a significant 
improvement over the SDSS PHOTO algorithm. However, our luminosity function of ETGs 
at the bright end is even shallower, particularly when we use the photometry to 1% of the 
sky brightness. 

Table 1 lists our measured luminosity densities for the Petrosian magnitudes and 
isophotal magnitudes to surface brightness of 25 mag and 1% of the sky background 
(denoted by 0p, and respectively) at several r-band luminosities. The bottom 
panel of Fig. [T3] shows the ratios of these galaxy luminosity densities to that measured by 
Blanton et al. (2003), denoted as 0p,Bianton, as a function of r— band luminosities. It is 
clear that the luminosity density ratios increase with the luminosity of ETGs. For ETGs 
with Mr = —23.5 mag, the ratios are 6.7 ± 0.46, 8.0 ± 0.40 and 10.5 ± 0.56, respectively; for 
ETGs with M, = -24 mag, they go up to 153 ± 14, 197 ± 16 and 259 ± 22, respectively. 
It demonstrates that the luminosity density of the brightest ETG calculated based on the 
SDSS catalogue has been seriously underestimated. Note that for the three magnitudes Mp, 
M25 and Mi%, the number of ETGs brighter than —23.5 mag is 347, 385 and 486, while 
the number of ETGs brighter than —24 mag is 41, 56 and 90, and so the number statistics 
are reasonably good. The underestimates in the luminosity density result in a significant 
underestimate of the integrated luminosity density. In particular, the integrated luminosity 
density down to Mr < —22.5 mag based on our measured Mp, M25, Mi% are about 20%, 
40%, 50% higher than that from the SDSS Petrosian luminosity, respectively. 

Next we estimate the stellar mass M* of the ETGs following Bernardi et al. (2010), 
utihzing the equation 

log M^/Mq = im7{g - r) - 0.406 - 0.4(Mr - 4.67) - 0.19;z, (3) 

where the g — r color is taken directly from the model color of SDSS photometrical catalog. 
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The equation has aheady taken into account the k-correction and evolution correction, and 
the initial mass function is assumed to be of the Chabrier (2003) form (see §2.4 of Bernardi 
et al. 2010 for more details). The stellar masses M^, p, M^, 25 and 1% are obtained using 
our measured luminosities Mp, M25, Mi%, respectively. 

Table 2 gives the stellar mass densities, 0m,, p, 0m,, 25 and 0m,, 1% at several stellar 
masses of ETGs for our measured magnitudes. The top panel of Fig. [15] shows the stellar 
mass function for massive ETGs. We also plot the stellar mass function from Blanton & 
Roweis (2007) and Bernardi et al. (2010) in the top panel of Fig. [15] for comparison. As 
can be seen, the slope of our measured stellar mass functions is shallower than those of 
Bernardi et al. (2010) and Blanton & Roweis (2007). Given that Bernardi et al. (2010) 
have already compared their stellar mass density with previous observational results, in the 
bottom panel of Fig. [15] we just compare our result with Bernardi et al. (2010). It shows 
that the stellar mass densities ratios of our measurements to that of Bernardi et al. (2010), 
written as 0m,, Bernardi? as a function of stellar mass. We can see from the bottom panel of 
Fig. [I5]that the ratios are larger for more massive ETGs. The ratios, 0m,, p /0m,, Bernardi, 
0Af„ 25 /0Af„ Bernardi and 0Af,, i%/0Af,, Bernardi are 1.2±0.06, 1.3±0.07 and 1.6±0.07 for ETGs 
with ~ 5 X lO^^M©; for ~ W'^Mq, the ratios are even larger, 2.1±0.40, 2.8±0.43 
and 4.2±0.58, respectively. Thus, at the high mass end the stellar mass densities have been 
underestimated by all previous works. 

We also plot the predicted stellar mass function by Guo et al. (2011) in the top panel 
of Fig. [IS Their result is based on the semi-analytic models of the galaxy population using 
the dark matter only Millennium Simulation. We can see that the high-mass tail of their 
stellar mass function slightly over- predicts the abundance in Bernardi et al. (2010), but 
still under-predicts our mass function. We return to this briefly in the next section. 
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5. SUMMARY 

In this work, we performed photometric analyses for a complete and homogeneous 
bright ETGs sample with 2949 early-type galaxies (M^ < —22.5 mag) in the redshift range 
of 0.05 to 0.15, taken from the catalog of SDSS DR7; all these galaxies have morphological 
classifications from the Galaxy Zoo 1 MGS. Based on our own sky background subtraction 
method, we measured the Petrosian and isophotal magnitudes, as well as the corresponding 
half-light radii. Comparing our measured luminosities and sizes to those from SDSS, we 
find that the SDSS pipeline significantly underestimates the luminosities and sizes for 
the brightest ETGs, leading to underestimates of the luminosity density and stellar mass 
density for bright ETGs. Our main results are summarized as follows. 

1. We find that our Petrosian magnitudes, and isophotal magnitudes to 25 mag/arcsec^ 
and 1% of the sky brightness are on average 0.16 mag, 0.2 mag, and 0.26 mag brighter 
than the SDSS Petrosian values, respectively. In the first case the underestimations 
are due to overestimations in the sky background by the SDSS PHOTO algorithm, 
while the latter two are also caused by an additional effect of reaching deeper 
photometry. Such underestimation is more severe as ETGs become more luminous. 
Our results also demonstrate that as we integrate to deeper surface brightness, we 
recover more luminosity of galaxies. 

2. We also find that the sizes of ETGs (half-hght radius r^o) measured by the SDSS 
algorithm are smaller than those measured by us. The largest r^o of bright ETGs in 
the SDSS catalogue is ~ 20 kpc and the mean is ~ 9 kpc, while our measured can 
be as large as 30 kpc and for a large fraction (~ 27%) of brightest ETGs, rso is larger 
than 15 kpc. An interesting result is that the underestimation of luminosities by SDSS 
is more strongly infiuenced by the size of ETGs than by their luminosities, indicating 
that for ETGs with the largest sizes, the luminosity underestimation by the SDSS 
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algorithm is more serious. In addition, we find that the slope in the size-luminosity 
relation at the bright end is much steeper than that found by Shen et al. (2003). 

3. Based on our selected complete and homogeneous sample of 2949 bright early-type 
galaxies, we construct the luminosity function. We find that the LF at the bright end 
is much shallower than those of Blanton et al. (2003) and Bernardi et al. (2010). The 
luminosity density at —23.5 mag (—24 mag) measured in this work is one order (two 
orders) of magnitude higher than that of Blanton et al. (2003). As a result, the ratios 
of the integrated luminosity density for bright galaxies (Mr < —22.5 mag) between 
those based on our measured Mp, M25, Mi% and the SDSS Petrosian luminosity are 
1.2, 1.4 and 1.5, respectively. Similarly, the stellar mass density of ETGs is a few 
tenths to a factor of few higher than that of Bernardi et al. (2010) for stellar mass 
from ~ 5 X lO^^M© to ~ IO^^Mq. Therefore, our method recovers substantially 
more luminosity for bright ETGs, which may alleviate the contradiction between 
hierarchical galaxy formation theories and current observations. 

Our results suggest very careful photometry needs to be performed to obtain the LF at the 
bright end which has significant impact on the stellar mass function for massive galaxies. 
Previous claims that the massive end of the stellar mass function has a weak evolution since 
redshift ~ 1 in comparison with the Galaxy And Mass Assembly (GAMA) survey (Baldry 
et al. 2012) may need to be re-examined. Furthermore, semi- analytical models may need to 
increase their coohng efficiency or decrease the AGN feedback efficiency for massive galaxies 
in order to match the shallower LF we found here. How does this affect the overall shape of 
the LF (especially around L*) is unclear and warrants further studies. 
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0.05 0.1 0.15 0.2 0.25 

redshift 

Fig. 1. — Three volume-limited subsamples selected from Galaxy Zoo 1 early- type galaxies 
in the redshift vs. absolute Petrosian magnitude plane in the r— band of SDSS DR7, shown 
in blue, cyan and magenta boxes, respectively. The total number of galaxies is 2949 ETGs. 
The red lines show the observed flux limits at r = 16 mag and 17.77 mag, respectively. 
The number of early-type galaxies in each subsample within the redshift range and lower 
luminosity limit are shown in the bottom right. 
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Fig. 2. — The sky coverage for our volume- limited bright early- type galaxies in the Aitoff 
projection in Galactic coordinates. The total sky coverage is about 9055 degree^. 
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Fig. 3. — Distributions of the SDSS Petrosian apparent magnitude in the r— band (left panel) 
and redshift (right panel) for 2949 early-type galaxies with Mr < —22.5 mag. The mean and 
median values for each distribution are shown in the top left. 
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Fig. 4. — Examples of objects that cannot be masked well by the automatic algorithm of 
SExtractor. The left panels show examples of color images of ETGs in a crowded field (top 
row), with extended stellar halo (middle row), and with nearby contaminated stars (bottom 
row). The middle and right columns show the corresponding masked images generated by 
SExtractor and by hand. The white areas indicate masked regions. 



- 30 - 




Fig. 5. — The difference Amp between the SDSS Petrosian magnitude and our measured 
Petrosian magnitude as a function of apparent Petrosian magnitude nip (top panel) and 
Petrosian absolute magnitude Mp (middle panel). The Spearman- Rank correlation coeffi- 
cients Ts are shown in the top left of the panels. The red data points with error bars are the 
median, lower (25 per cent) and upper (75 per cent) quartiles for binned galaxies. The bin 
width is 0.4 mag for the apparent Petrosian magnitude except the last bin which has a width 
of 0.6 mag to include all the remaining objects (top panel); the bin width is 0.2 mag for 
the absolute magnitude except the last bin which has a width of 0.3 mag to include all the 
remaining objects (middle panel). The bottom panel shows the histogram of the magnitude 
difference with the mean and median values indicated in the top right. 
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Fig. 6.— The distributions of the difference between the SDSS DR8 and SDSS DR7 IPet- 

rosian magnitudes (left panel), and the difference between our measured Petrosian magnitude 
and the SDSS DR8 (right panel). The mean and median values for each distribution are 
shown in the top left of the panel. 



- 32 - 




Fig. 7. — The difference Am25 between the SDSS Petrosian magnitude and our isophotal 
magnitude to 25 mag/arcsec^ as a function of apparent 17125 (top panel) and absolute M25 
(middle panel) isophotal magnitudes. The Spearman-Rank correlation coefficients are 
shown in the top left of the panels. The bottom panel shows the histogram of magnitude 
difference distribution with the mean and median values indicated in the top right. See 
Fig. |5]for an explanation of the red data points with error bars in the top and middle panels. 
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Fig. 8. — The difference Ami% between the SDSS Petrosian magnitude and our isophotal 
magnitude with surface brightness hmit at 1% of sky brightness as a function of the apparent 
mi% (top panel) and absolute Mi% (middle panel) magnitudes. The Spearman-Rank corre- 
lation coefficients Ts are shown in the top left of the panels. The bottom panel shows the 
histogram of magnitude difference distribution with the mean and median values indicated 
in the top right. See Fig. [5] for an explanation of the red data points with error bars in the 
top and middle panels. 
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Fig. 9. — The histograms of half-hght radii (^50) measured by different photometric methods. 
The black and green dotted hnes represent the Petrosian r^Q distributions based on SDSS and 
our own measurements, respectively. While the red solid and blue dot-dashed lines show the 
rso distributions based on isophotal measurement to 25 mag/arcsec^ and 1% of sky bright- 
ness, respectively. The median and mean r^Q values (in units of kpc) for each distribution are 
shown at the top right of the figure. Log-normal (see eq. [2]) fits are shown for each distribu- 
tion with corresponding colors; the medians and dispersions are f = 8.58, 10.28, 10.80, 11.82 
kpc and ainr = 0.24,0.32,0.31,0.32 for r5o,sdss7, ^50,?, r5o,25, ^50,1% respectively. 
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Fig. 10. — The magnitude differences between the SDSS Petrosian and our measured Pet- 
rosian, isophotal magnitudes with surface brightness hmits at 25 mag/arcsec^ and 1% of sky 
brightness respectively, as a function of our measured corresponding half-hght radius in log- 
arithm. The red data points with error bars are the median, lower (25 per cent) and upper 
(75 per cent) quartiles for galaxies in bins of width 4 kpc in half-light radius except the first 
bin which has a width of 3 kpc and the last bin which has a width of 10 kpc to include all 
the remaining objects. The mean, median magnitude differences and the Spearman-Rank 
correlation coefficient are shown in the top left of each panel. 
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Fig. 11. — The upper panels show four (typical) examples of bright ETGs located in crowded 
fields, and the lower panels are examples of those contaminated by bright nearby stars. 
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Fig. 12. — Examples of bright ETGs with extended halos and corresponding surface bright- 
ness profiles. The red solid line shows the best fit Sersic law convolved with the point spread 
function. The dotted and dashed lines in the middle panels represent the surface brightness 
of 25 mag/arcsec^ and 1% of sky brightness respectively. The corresponding Sersic law in- 
dex n is shown in the top right of each panel. The bottom panels show the residuals of the 
observed surface brightness from the Sersic model. 
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Fig. 13. — The grayscale representation of the size-luminosity relation for 2949 bright ETGs. 
The four absolute magnitudes used are the SDSS Petrosian magnitude (top left), our own 
Petrosian magnitude (top right), and aperture magnitudes to 25 mag/arcsec^ (bottom left) 
and 1% of sky brightness (bottom right), respectively. The red line in each panel is the 
best power-law fit for the size-luminosity relation. The power-law index a [t^q oc L") is 
shown in the bottom right of each panel. The rightmost vertical grayscale bar refiects the 
corresponding number of galaxies. 



- 39 - 




Fig. 14. — The top panel shows the r— band luminosity function at the bright end calculated 
based on our measured luminosities for 2949 ETGs brighter than Mr < —22.5 mag. The 
green triangles, red circles, and blue solid squares represent the LF calculated by using our 
measured Petrosian magnitude, isophotal magnitudes with surface brightness limits at 25 
mag/arcsec^ and 1% of sky brightness respectively. The magenta dot-dashed and black solid 
lines give the fits by Blanton et al. (2003) and Bernardi et al. (2010). The points for objects 
brighter than the luminosity shown by the cyan dotted line are not affected by photometry. 
The bottom panel shows the galaxy luminosity density ratios as a function of the r— band 
luminosity. The green triangles, red circles, and blue solid squares show the ratios for our 
measured Petrosian, isophotal luminosities to 25 mag/arcsec^ and 1% of sky brightness to 
that from Blanton et al. (2003). 
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Fig. 15. — The top panel shows the stellar mass function for massive ETGs. The stellar 
masses are estimated from our measured Petrosian (green triangles), isophotal luminosities to 
25 mag/arcsec^ (red circles) and 1% of sky brightness (blue solid squares), respectively. The 
magenta dot-dashed and black solid lines represent the best-fit stellar mass functions from 
Blanton & Roweis (2007) and Bernardi et al. (2010), respectively. The navy blue dashed 
line represents the prediction from the semi-analytic models of the Millennium Simulation 
of Guo et al. (2011). The points for objects more massive than the stellar mass shown by 
the cyan dotted line are not affected by the photometry. The bottom panel shows the ratios 
between our measured stellar mass densities and those of Bernardi et al. (2010) as a function 
of stellar mass; the symbols are the same as in the top panel. 



Table 1. Luminosity densities for the Petrosian magnitude (0p) and isophotal magnitudes 
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,551 


634.400±31, 


,718 


645.401±31, 


,992 


-22.95 


443.715±26. 


,526 


522.556±28. 


,787 


573.894±30. 


,168 


-23.05 


335.537±23. 


067 


425.38a±25, 


,973 


445.300±26, 


251 


-23.15 


244.257±13. 


,778 


286.031±21, 


,298 


340.198±24, 


,652 


-23.25 


189.786±12. 


,098 


218.191±18. 


,601 


262.195±20. 


,391 


-23.35 


127.568±9. 


,925 


159.517±15. 


,905 


203.522±17. 


,965 


-23.45 


92.339±8. 


,424 


118.44S±9. 


,555 


146.683±15, 


,252 


-23.65 


73.341±10. 


785 


82.101±7, 


963 


114.717±13, 


132 


-23.65 


49.505±8. 


,860 


64.004±7, 


,017 


99.011±12, 


.530 


-23.75 


35.526±5. 


226 


47.672±8. 


,695 


66.697±8, 


.242 


-23.85 


20.151±2, 


,954 


24.958±6, 


,438 


33.346±8, 


,371 


-23.95 


18.670±2. 


,945 


22.846±3. 


195 


31.383±5, 


.785 


-24.05 


12.094±2, 


,348 


16.537±2. 


,857 


20.738±3, 


.073 


-24.15 


4.397±1. 


,394 


7.915±3, 


.681 


17.561±3, 


.025 


-24.25 


3.582±1. 


,923 


4.075±1, 


.923 


4.776±2, 


.221 


-24.35 


1.194±1. 


111 


2.262±1, 


.071 


4.234±2, 


.784 



^The luminosity density at Mr = —22.75 mag is affected by the photometry. 



Table 2. Stellar mass densities for the Petrosian magnitude (0M«,p) and isophotal 
magnitudes to 25 mag/arcsec^ (0m., 25) and 1% of the sky background (0m., i%)- 



log M,/M0 


(10-'^ Mpc-^ dex 




"^M, , 25 
(10-'^ Mpc-^dex 




'I'M, , 1% 
(10-^ Mpc-^dex 




11.425=' 


1503.490±69. 


055 


1481.4!)()±68, 


,547 


1437, 


480±(i7, 


522 


11.475 


1078.120±58. 


476 


1239.460±62. 


,699 


1364, 


,140±65, 


,777 


11.525 


876.429±52. 


,723 


964.438±55. 


,307 


982. 


,773±55. 


,830 


11.575 


685.742±46. 


636 


722.413±47. 


,867 


880, 


,096±52, 


,833 


11.625 


370.374±34. 


,274 


476.719±38. 


,884 


568. 


,396±42. 


,459 


11.675 


267.674±28. 


369 


333.723±31. 


,676 


421. 


,713±36. 


,572 


11.725 


202.038±24. 


,976 


201.688±25. 


,292 


260. 


,721±27. 


,998 


11.775 


128.623±19. 


,665 


161.351±22. 


,622 


192. 


,372±24. 


,181 


11.825 


95.344±17. 


390 


113.685±18, 


,989 


135, 


768±20, 


,758 


11.875 


55.620±12. 


,932 


62.573±13, 


.716 


80, 


,665±13. 


,755 


11.925 


36.619±9. 


,903 


38.211±8, 


.884 


56. 


,210±13. 


,499 


11.975 


19.166±7. 


968 


26.762±8, 


.448 


39. 


lOOill. 


,500 


12.025 


11.204±6. 


141 


16.896±7, 


.393 


23.a64±8, 


.333 


12.075 


5.627±4. 


,233 


8.355±5, 


.495 


12.6a4±6, 


.768 



^The stellar mass density at log M* /Mq = 11.425 is affected by the photometry. 



